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Abstract 

This is not a new result. Purpose of this work is to describe a method to search 
the analytical expression of the general real solution of the two-dimensional 
Laplace differential equation. This thing is not easy to find in scientific liter- 
ature and, if present, often it is justified with the assertion that an arbitrary 
analytic complex function is a solution of Laplace equation, so introducing the 
condition of complex-differentiability which is not really necessary for the exis- 
tence of a real solution. The question of the knowledge of real exact solutions 
to Laplace equation is of great importance in science and engineering. 
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Consider the 2D Laplace equation dxxU + dyyU — for a function u = u{x, y) 
defined on a region 17 C M^. If one, e.g. for physical or engineering reasons, 
is interested in the symbolic (analytical) expression of the general real solution, 
usual technical good literature (see [4] for a general reference) gives three cate- 
gories of answers: 

a. theorems of existence and uniqueness for boundary value problems, very im- 
portant results but often without a symbolic expression of the solution; 

b. applications of the method of separable variables, which often gives a sym- 
bolic solution of the differential equation but only in the restricted range of 
functions of the form F{x)G{y); 

c. an information, unfortunately often not completely justified, about the fact 
that an analytic (in the sense holomorphic) complex function F{z) is harmonic. 
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that is AF{z) = 0, where is the origin of the complex plane. 



In particular, in the latter case the function F must be C-derivable, that is 
dxF = —idyF, which is a condition very stronger than the differentiability of the 
real and imaginary part; therefore, complex analytic functions are a restricted 
range of all the possible solutions of the Laplace equation. 
There is another way to search the form of the solution: the use of a symbolic 
mathematical software. If we ask Mathematica to solve dxxU + dyyU = 0, the 
built-in function DSolve answers 



where = — 1 and F, G are arbitrary functions. The relative tutorial gives 
no more informations about the nature of these arbitrary functions. Note that, 
in general, F(y + jx) and G{y — jx) are complex numbers. Previous formula ([T]) 
is present in ([6]), but only as exercise to show that it is a solution of Laplace 
equation. 

In[l];= LselaceEaiLatlon = D[u[x, y] , {x, 2 }] + D [u[x, y] , {y, 2}] == 0; ^ 
In[2]i= DSolire [LaplaceEguation, ti[x, y] , {x, y}] ^ 

Out[21= {{u[x, y] ^ C[ll [±x + Yl + C[21 [-ix + Yl }} 3. 

This general solution contains two arbitrary functions, c [ 1] and c [2 ] . 



Hence, what is the mathematical form of the real general solution of the Laplace 
equation? In this work I try to find it with a method which is not standard in 
the usual literature. 

Consider the laplacian operator A = dxx + dyy as an algebraic object defined 
on the set of twice partially differentiable real functions defined on R^. We can 
define the object [dxf = dxx, that is {dxf f = dxdxf for every / G C'^{R^,R). 
Then it can be easily stated the following identity: 



because on C^(R'^,R) the two single elementary operators dx and dy commute: 



{dx){jdy)f = {jdy){dx)f. Define Dj = {dx + jdy) and = {dx - jdy). Two 



similar operators are defined, in the discussion of the Cauchy-Riemann equa- 
tions, in the first edition of ( 5 ). 

Let be M e C^(K.^, M); then, being D-jDjU = DjD^jU, u is a solution of Au = 
if DjU = or D-jU — 0. Consider the first case: 



u{x, y) = F{y + jx) + G{y - jx) 



(1) 



Figure 1: Output from DSolve. 



{dx+ jdy){dx - jdy) 



(2) 



dxU + jdyU = 



(3) 
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It is easy to apply the method of characteristics (see [5]) for its symboHc reso- 
lution. The ordinary differential equation resolved by the characteristics is 



y'(^)=j (4) 

which is satisfied by the family y = jx + c, with c arbitrary (complex) constant. 
Along each characteristics the equation ^ becomes (note that the complex 
function jx + c is R-derivable with respect to the real variable x) 

dxiu{xjx + c)) ^0 (5) 

so for every choice of c there is an arbitrary (complex) constant G, depending on 
c, such that u{x, jx+c) = G(c) . But for every value of y there is a characteristics, 
namely the line corresponding to the constant c = y — jx, passing for the point 
{0,y) of K^), so if we take a (complex) function G whose real and imaginary 
part are twice M-differentiable on R^, the general real solution of ^ is 

ui = Re [G{y - jx)] (6) 

In the same manner, one finds that U2 ~ F{y + jx) is the general complex 
solution of the second case, D^ju = 0. 

Now, suppose Au = DjD-jU — 0. If D^ju = or DjU — 0, the case has already 
been treated. Then, suppose e.g. D_jU ^ 0. Therefore Dj {D_ju) = 0, so from 
^ there is a function G, with real and imaginary part twice R-differentiable on 
R^, such that D^ju = G{y — jx). Now, to find u, we have to resolve the non 
homogeneous PDE 

dxU - jdyu = G{y - jx) (7) 

But the homogeneous equation D^jU — has general solution Uh = F{y + jx). 
We try to find a particular solution Up oi Let be Up{x, y) = a{x, y)+jb{x, y). 
Then 

D^jUp = {dxtt + dyb) + j(-dya + d^b) (8) 
and the two following first order pdes must hold (Gr — M.e[G], Gj — Im[G]): 

d^a + dyb = Gr, -dyU + d^b = Gi (9) 
A possible way to find a solution is the choice 



d^a = ^Gr, -dya = ^Gi 
^Gi, dyb = i( 



dxb = -Gi, dyb = -Gi 



which can be resolved by quadrature. Let be A the complex function so defined: 
A{z) = A{x + jy) = a{-y, x) + ]b(-y, x) (10) 
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Note that A{y — jx) = a{x, y) + jb{x, y), so the function Up{x, y) = A(j/ — jx) 
is a particular solution of ([7]). We have shown that if Au = 0, then u has the 
form ll]). 

In conclusion, the general real solution of the Laplace equation is 

u{x, y) = Re [F{y + jx) + G{y - jx)] (11) 

Note that for our considerations it is not necessary the C-differentiability of F 
and G. Therefore, the assertion about Laplace equation, made by Mathematica 
and by some texts on complex analysis, is now more clear when transposed in 
the context of real functions: 

the general real solution of Laplace equation dxxU + dyyU ^ Q is u{x,y) = 
Ee [F{y + jx) + G(y — jx)] where F and G are arbitrary complex functions such 
that the real part of F(y + jx) + G{y — jx), considered as function of {x, y), is 
twice R-dijJerentiable. 

Example. Suppose we want to find a solution of the differential problem dxx"^ + 
dyy'i' = on the domain = {{x,y) x > 0,0 < y < h}, where h > 0, with 
boundary condition 5* = on dfl. 

Such a problem arises in fluid dynamics when 'J is the stream function of a flow 
in n, semi-infinite bay, with velocity field {dy'i', —dx'i') (see [5]). 
Note that the problem has not a unique solution. The trivial function 5* = 
is solution, and for every 5" solution of the problem, fc^* is solution too for 
every constant k. We try to construct a not trivial solution from the general 
formula (fTTj) . Being '^{0,y) — for every y G [0,h], a possible solution could 
be \1/ = sin(x)/(y) with / such that /(O) = f{h) = 0. One knows that for the 
complex function cos(ja; + y) the following identity holds: 

cos(ja; + y) ~ cosh(x) cos(y) — j sinh(a;) sin(y) (12) 

The function j cos{jx + y) is of type F{y + jx), and its real part is a twice 
R-differentiable function. For the boundary conditions, note that it is sufficient 
to apply a simple change of variable jx + y if (-^^ f)' "^^icre n G N. The 
real part becomes 

vl.(x,.)^sinh(l^)sin(!^) (13) 

that is a real solution of the boundary problem. The same problem, with the 
same solution, is treated in ([3]), where separation of variables is applied. 

Application. Find the general real solution of the PDE 

dx.U ~ dyyU = (14) 

This equation appears in the description of the shear stress of a plane flow 
with velocity field (dyU, —dxU). We can repeat previous argument with dxx — 
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dyy = [dx + dy){dx ^ Qy) . But It Is 111016 liiterestiiig if we apply the variables 
transformation X ~ x, Y = jy, where = —1. Then (|14p becomes dxxU + 
dyyU = where U = U{x{X),y{Y)), which is the previous Laplace equation. 
Its general real solution is U{x,y) = Re [F{Y{y) +jX{x)) + G{Y{y) - jX{x))], 
therefore, being Y + jX = jiy + x) and Y — jX — j{y — x), we can write 

U{x,y)=Re[F{y + x)+G{y~x)] (15) 
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